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Abstract 

We have carried out an inclusive measurement of production in two- 

photon processes at TRISTAN. The mean -y/s was 58 GeV and the integrated 
luminosity was 199 pb"-*^. High-statistics Kg samples were obtained under such 
conditions as no-, anti-electron, and remnant-jet tags. The remnant-jet tag, in par- 
ticular, allowed us, for the first time, to measure the cross sections separately for 
the resolved-photon and direct processes. 
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1 Introduction 



Excess over theoretical predictions was reported in references [||, H, concerning charm pair 
production in two-photon processes. These references extensively discussed possibilities to 
explain the excess by increasing the predicted cross sections of two-photon processes. In 
order to help sorting out these possibilities, we have carried out an inclusive measurement 
of K^{K^) mesons in two-photon processes. With the point-like process (direct process 
0]) where e'^e^ e'^e~ss is strongly suppressed due to the small value of {Qs is s 
quark's charge), Kg^s, particularly in the high-P^ region, come mainly from cc production 
The use of Kg^s also enhances the sensitivity to two-photon processes at low Pt- 
There are additional advantages in this Kg analysis, which are high statistics and 
low systematics, thanks to high acceptance and use of only a central tracking device, 
respectively. Moreover, using the low angle calorimeter [Q, we could obtain the cross 
section for the resolved photon process separately. We can thus test some of the 
possibilities suggested in reference 0: we proposed to use a low charm quark mass (1.3 
GeV), the next-to-leading-order correction (NLO), and the intrinsic parton Pt inside 
photon in order to explain an excess observed in the cc cross section especially in the 
high-Pr region. 

Notice also that this analysis is hardly affected by t (superpartner of the top quark) 
pair production discussed in references 0], even if present, since the high-Q^ decay 
of c — s> K'^X smears the Pt distribution of the Kg's and thus diminishes the sensitivity 
to the t pair production. 
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2 Event selection 



The data used in this analysis were obtained with the TOPAZ detector at the TRISTAN 
e^e~ colhder, KEK||, The mean a/s was 58 GeV and the integrated luminosity was 
199 pb^^. A forward calorimeter (FCL) , which covered 0.98 < |cos^^| < 0.998 {9 is the 
polar angle, i.e., the angle with respect to the electron beam), was installed in the course 
of the experiment. The FCL was made of bismuth germanate crystals (BGO), and was 
used to anti-tag the beam electrons (positrons) and to tag hadrons (remnant-jets) 
for this study. We could thus select collisions of almost-real photons including resolved 
photon processes 0. The integrated luminosity of the data with the FCL detector was 
175 pb-^ 



A description of our trigger system can be found in reference [|T0|. The requirement 
for the charged track trigger was two or more tracks with an opening angle of greater than 
45-90 degrees. The Pt threshold for charged particles was 0.3-0.7 GeV, varied depending 
on beam conditions. 

The event selection criteria were as follows: three or more charged particles {Pt > 0.15 
GeV, I cos 6*1 <0.77), the invariant mass {Wvis) of visible particles (| cos^^| <0.77) had to 
be greater than 2 GeV, the event- vertex position had to be consistent with the interaction 
point, and the visible energy had to be less than 25 GeV. In total, 220378 events were 
selected. 

3 Monte-Carlo simulation 

In order to estimate the acceptances and backgrounds in this analysis, we used the fol- 
lowing Monte-Carlo simulation programs. For the generation of single-photon-exchange 



hadronic events, we used JETSET6.3|llTl] with the parameter values given in reference 
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|T^. Details concerning the event generation of direct as well as resolved-photon and 
vector meson dominance ( VDM) processes can be found in references |1|, 0, |T^ . Here, we 
just note the following points. For cc generation, we used the current charm quark mass 
of 1.3 GeV to calculate cross sections for point-like processes and a constituent charm 
quark mass of 1.6 GeV for hadronization procedure, and made the next-to-leading order 
(NLO) correction, whose details can be found in references |l|, ^, |l^. Light-quark genera- 
tion was carried out by using the lowest order (LO) formula with a P™" cut of 2.5 GeV. 
We used the parton density functions by Levy-Abramowicz-Charchula set-1 (LACl) [|l^ 
for the resolved-photon process. Generated events were processed through the standard 



TOPAZ detector simulation program [0, in which hadron showers were simulated with 
an extended version of GHEISHA 7 |jl3[ . Its data on hadron interactions with nuclei had 
been updated to fit various experimental cross sections. 

Using the above-mentioned Monte-Carlo simulations, the trigger efficiency for the 
sum of the direct and resolved photon processes was estimated to be 79%, 97% of which 
represented charged trigger events. The event-selection efficiency after the trigger was 
obtained to be 80%. 

4 Tagging conditions 

The tagging conditions were as follows. For anti-electron tagging, there had to be no 
energy deposit of more than OAEb in | cos^^| < 0.998 (anti-electron tag or anti-tag), where 
Ef, is the beam energy. This selected events from collisions of almost real photons, for 
which the accuracy of the equivalent photon approximation was expected to be reliable 
at the 1% level. When the energy cut was lowered, mis-anti-tag due to beam remnant 
hadrons (remnant-jets) became significant, as predicted by the Monte-Carlo simulations. 
The energy distribution of the maximum-energy cluster in FCL is shown in Figure 



The horizontal scale is normalized at the beam energy. The Monte-Carlo predictions for 
single-photon-exchange, VDM, and resolved-photon processes are shown by histograms. 
The peak around 1 was explained by the direct processes. This implies a possibility to tag 
the resolved photon process by requiring, for instance, 500 MeV < EpcL < 0.25i?b, where 
EpcL is the energy deposit in the FCL (remnant-jet tag or rem-tag). We did not use any 
hadron shape information because of large segmentation and lack of tracking information. 
The yield of the remnant-jet tag events agreed with our Monte-Carlo simulation within 
5% level. Our analysis is the first trial that uses this tag. An event selection without these 
two tags is called, hereafter, "no-electron tag" or "no-tag" . The fractions of electron and 
remnant-jet tag events to no-tag ones were obtained to be 2.4 and 47% of the selected 
events, respectively. 

In the Monte-Carlo simulations, we used the equivalent photon approximation with 
the photon flux expression given in reference [0. We set the Q"^ limit at the smaller of 
Pj-^q + ^l and the anti-tag limit [2E'^{l — x^){l — cos9max) '■ = 0.4, = 3.2 degrees], 
where Pt^q and mg are the transverse momentum and the mass of a quark, respectively 

The tagging efficiency of the remnant-jet tag for the resolved photon process was 
estimated to be 72% without assuming FCL noise which will be described later. We tried 
two ways of generating remnant partons: one along the beam direction, and the other 
using a Gaussian distribution of P-r-width 0.44 GeV with respect to the beam axis. Two 
methods differed in acceptance only by 3%. On the other hand, the tagging efficiency of 
the remnant-jet tag for the direct process was estimated to be 10.8% without assuming 
FCL noise. 
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5 Analysis 

The charged-track selection criteria for the analysis were as follows: for each track 
had to be greater than 0.15 GeV, dE/dx had to be consistent with the pion assumption 
(x^± < 10), I cos^l had to be less than 0.77, and the closest approach to the interaction 
point in the XY-plane (perpendicular to the beam axis) had to be greater than 1 cm. Using 
these selected tracks, we looked for opposite-sign pairs with an opening angle less than 90 
degrees, and carried out secondary vertex reconstructions three-dimensionally. Finally, we 
demanded these pairs to be consistent with the assumption that they came from the event 
vertices with a flight length larger than 3 cm. The invariant-mass distributions of these 
candidate tt^tt' pairs are plotted in Figure]^ for the three tagging conditions, respectively. 
These invariant-mass spectra were fitted with the sum of a second-order polynomial and 
a Gaussian distribution and the peak entries were obtained for no-, remnant-jet, and 
electron tags to be 893 ± 34, 364 ± 22, and 75.8 ± 9.5 -R's's, respectively. The peak 
position and the width were consistent with the detector simulation. In order to derive 
the differential Pt cross sections, we divided Pt into ten bins, as shown in Table ^ Notice 
that even the lowest-statistics bin gave a 5.1-cr Kg peak (in the no-electron tag case). 

6 Background subtractions 

Single-photon-exchange process was a largest background especially for high-P^ K^. This 
can be reduced when a cut was applied on the total visible energy. We, however, did not 
carry out this, because we did not want to reduce the acceptance for high-Py Kg. The 
contamination from the single-photon-exchange process was estimated and subtracted 
using the Monte-Carlo simulation, on a bin by bin basis. The background fractions for 
no-tag were 6.0, 6.1, 6.5, 8.1, 10.0, 11.0, 14.3, 17.3, 43.4, and 48.1%, respectively, for the 
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Pt bins shown in Table [1|, being strongly Pt dependent. We estimated these fractions for 
anti- and remnant-jet tags. They were consistent with the above values within statistical 
errors. The background from beam gas interactions was estimated using the off-vertex 
events in the beam direction: there was a vacuum leak in the beam pipe for some period. 
The beam gas contribution for no-tag was 7.8% on the average, and was subtracted from 
the data. For anti- and remnant-jet tags, the above value became slightly large (9.8%), 
i.e., electron-tag sample was free from the beam-gas background. FCL noise hits were 
studied by analysing random trigger and Bhabha events. The probability of noise hits with 
EpcL >0.5 GeV was estimated to be 12.1%, while for hits with EpcL > 0.4i?;,, it reduces to 
0.1%. The FCL noise was also related to the vacuum leak. In the Monte-Carlo simulations, 
we added noise hits randomly in accordance with the observed noise hit probability in 
order to reliably estimate the tagging efficiencies. Using this, the tagging efficiencies of 
the remnant-jet tag for the resolved-photon and direct processes were estimated to be 
75% and 21%, respectively. 

7 Systematic errors 

The systematic errors were estimated, bin by bin, as follows. For the trigger, we added 
some extra noise hits in the tracking chambers in the simulations. For the event selection 
and the Kg reconstruction, we changed the cut values and evaluated systematic errors 
as the cross-section differences. We also changed the pulse-height threshold in the TPC 
simulation to evaluate the effects on its tracking efficiency. We added the obtained sys- 
tematic errors quadratically, on a bin by bin basis. The average over the Pt bins of the 
systematic errors was 12%, of which the cut dependence in the event selection was the 
dominant source. These systematic errors were quadratically added to statistical errors. 
We also checked the acceptance ambiguity due to the parametrization dependence of 
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the resolved-photon processes, by comparing the LACl |]T6| and Drees-Grassie [DG] [|T^ 
parametrizations. The acceptance difference was estimated to be 5.9%, which is small 
compared to the systematic errors shown above. 



8 Results 

The Pt differential cross sections were obtained from the number of reconstructed Ks^s 
in each bin and its corresponding efficiency estimated with the Monte-Carlo simulations 
described previously. They are listed in Table and plotted in Figures ^ (a) - (e) for the 
three tagging conditions and two subtraction schemes, respectively. Figure]^ (a) is for anti- 
electron tag events. In the remnant-jet tag events, the Monte-Carlo simulation predicted a 
significant contamination from VDM events. The tagging efficiency for the VDM process 
was estimated to be ~61%, slightly smaller than that of the resolved-photon process. In 
addition there was a large ambiguity in the cross section of the VDM process. Therefore 
we calculated the cross sections using two subtraction schemes. Figures ^ (b) and (e) were 
obtained by subtracting the VDM contribution predicted by the Monte-Carlo simulation 
(VDM subtraction) for remnant-jet and anti-remnant-jet tags, respectively. Here, the 
"anti-remnant-jet tag" cross section was obtained by subtracting the remnant-jet tag cross 
section from that of the anti-electron tag. Figures |^ (c) and (d) were obtained without 
the VDM subtraction. In Figures ^ (b) and (d), the contribution of the direct process was 
subtracted, since the uncertainty in the prediction of this process was considered to be 
small. The histograms in Figures |^ (a) - (e) are the Monte-Carlo predictions: the cross- 
hatched, singly-hatched, and open areas are predictions for the direct, resolved-photon 
(LACl), and VDM processes. 
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9 Discussions 



The fraction of charm events was studied using the above-mentioned Monte-Carlo simu- 
lations. We found 55% of these events with Pt{Ks) > 1.4 GeV were of charm origin. On 
the other hand, only 30% of the events with charged tracks of > 1.4 GeV were from 
cc pairs. Also, the Monte-Carlo simulations predicted that 70% of these high-Py charm 
events originated from the direct process. In this study, we derived six types of cross 
sections using different tagging conditions and subtraction schemes. We can therefore 
separately compare each cross section with the theoretical prediction for each process. 

Firstly, about 30% of the high-P^ events can be explained as the electron-tagged 
events (see Table |1|). Secondly, for the anti-tag cross section [Figure (a)], the agreement, 
especially in the lower- Pj- region, is reasonably good, considering the ambiguities due to 
the VDM process. They are, however, higher than the theoretical predictions in the 
high-Pr region (Pr > 1.2 GeV) hj 2.2a. 

The cross section with the anti-remnant-jet tag is consistent with the predictions of 
the direct and VDM processes [Figures^ (c) and (e)]. The hypothetical t pair production 
with mi=15 and m;y=12.7 GeV [jl], is expected to increase these cross sections by 1.7, 
1.1, and 0.3 pb/GeV in the highest three Pr bins, respectively, which are smaller than 
the experimental errors. Therefore we could not discuss this hypothesis by this analysis 
result. Note also that this hypothesis turned out to be inconsistent with the recent search 
by the VENUS collaboration 0. 

The dominant source of the discrepancy lies in Figures |^ (b) and (d), i.e., in the 
remnant-jet tag sample (we considered that this sample was dominated by the resolved- 
photon and/or VDM processes): the spectrum looks harder than the prediction of the 
resolved photon process. The histograms in Figures |^ (b) and (d) have already been 
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corrected by the NLO factor (0.5Pt,c+0.54; Pt,c is a Pt of a charm quark) [|T], 0, p!5| . 
This factor, which is due to the presence of off-shell gluon emissions from resolved quarks, 
is large especially in the high Pt region. The measured spectrum is harder than the 
LO prediction, and is close to the spectrum of the direct process. Our data suggest the 
importance of the NLO correction. The NLO correction to the light quark events, which 
is absent from our present Monte-Carlo generator, is also considered to be necessary. In 
addition, the intrinsic parton Pt (we used 0.44 GeV with a Gaussian distribution) is 
necessary. In summary, the theoretical prediction agrees with our data to a reasonably 
good accuracy, justifying our parametrization of the Monte-Carlo generation. 

In order to check if the parton density functions have anything to do with the discrep- 
ancy, we compared our remnant-jet-tag data [Figure |^ (d)] with the predictions from six 
sets of parametrizations by Hagiwara, Tanaka, Watanabe, and Izubuchi [WHITl-6] p!8| . 
A systematic analysis on gluon distributions can be carried out using these parametriza- 
tions. The results are shown in Figures ^ (a) and (b), where the histograms are the 
predictions by the WHIT 1-6 parametrizations with P^™'s of 2.0 and 2.5 GeV. Notice 
that the predictions are rather sensitive to the P™" cut and the VDM contribution, es- 
pecially in the low-P^- region; P^ < 1 GeV. There are some possible combinations which 
reproduce the high-Pr cross sections well. It is, however, necessary to improve the pre- 
diction of the VDM process. 

10 Conclusion 

We carried out an inclusive measurement of K^{K^) in two-photon processes at TRISTAN. 
The mean ^/s was 58 GeV and the integrated luminosity was 199 pb^^. High- statistics 
Kg samples were obtained under such conditions as no-, anti-electron, and remnant-jet 
tags. Especially with the remnant-jet tagging, we could unambiguously extract the con- 
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tribution from the resolved photon process. Comparisons with the theoretical predictions 
were carried out. Our results agreed with the theoretical predictions with a low charm 
mass (mc=1.3 GeV), intrinsic parton Pt inside photon (crpy=0.44 GeV), and the NLO 
corrections, i.e., those obtained by the previous results of the D*^ analysis [H, |]. 
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Table |t|, R. Enomoto et al., Physics Letters B. 



ta£f cond 


no-taff 


anti-ta£f 


rem-tag 


anti-rem-tag 


rem-ta£f 

X Vvxxx u caj 


aiiti-rem-tcLg 


V ly L\L o Li U L . 






yes 


no 




yes 


Pt range 






cross sections 






(GeV) 






(pb/GeV) 






0.3-0.55 


889± 248 


810± 247 


43± 52 


767± 247 


613± 237 


197± 282 


0.55-0.65 


839± 166 


783± 168 


53± 42 


731± 170 


520± 161 


263± 202 


0.65-0.75 


703± 114 


593± 109 


99± 34 


494± 108 


556± 117 


37± 137 


0.75-0.9 


415± 73 


371± 70 


34± 20 


337± 69 


241± 62 


129± 74 


0.9-1.05 


236± 32 


213± 32 


41± 14 


172± 33 


161± 33 


52± 42 


1.05-1.2 


162± 24 


116± 22 


14± 10 


102± 23 


51± 20 


65± 28 


1.2-1.4 


70± 12 


62± 12 


16± 7 


46± 14 


30± 10 


32± 15 


1.4-1.7 


37. 7± 8.3 


30.6± 7.6 


15.8± 5.3 


14.8± 8.0 


19.1± 6.1 


11. 5± 8.3 


1.7-2.5 


7.2± 2.2 


4.8± 2.0 


4.7± 1.8 


0.1± 2.5 


4.8± 1.9 


0.0± 2.5 


2.5-5 


1.78±0.73 


1.63±0.70 


0.98±0.64 


0.65±0.93 


0.98±0.64 


0.65±0.93 



Table 1: Differential cross section of K^lK'^) versus Pt (GeV), da/dPx (pb/GeV), for 
I cos^l < 0.77. Six cases are listed: no-tag , anti-electron tag, remnant-jet tag with VDM 
subtraction, anti-remnant-jet tag without VDM subtraction, remnant-jet tag without 
VDM subtraction, and anti- remnant-jet tag with VDM subtraction, which are described 
in the text. 
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Figure captions 



Figure 1: Distribution of the energy fractions (normalized at the beam energy) of the 
maximum-energy clusters in FCL. The points with error bars are experimental data. 

The histograms are the prediction by the Monte-Carlo simulation; the cross-hatched area 
is single-photon-exchange process, the singly-hatched one is VDM, and the open one is 
resolved-photon process. 

Figure 2: Invariant-mass spectra of tt+tt^'s. The solid histogram is for the no-electron 
tag, the dashed one is for the remnant-jet tag, and the dotted one is for the beam-electron 
tag. 

Figure 3: Differential cross section of K^{K'^) versus Pt (GeV), da/dPr (pb/GeV), 
for |cos^^| < 0.77. Five cases are plotted: (a) anti-electron tag, (b) remnant-jet tag 
with the VDM and direct process subtraction, (c) anti- remnant-jet tag without the VDM 
subtraction, (d) remnant-jet tag without the VDM subtraction, and (e) anti-remnant- jet 
tag with the VDM and direct process subtraction, as described in the text. Processes 
which we expected to show are; (b) resolved-photon, (c) direct and VDM, (d) VDM and 
resolved, and (c) direct processes. The histograms arc the theoretical predictions which 
are described in the text. The open area is for the VDM, the singly-hatched one is for 
the resolved photon process, and the cross-hatched one is for the direct process. 

Figure 4: Differential cross section of K^{K^) versus Pt (GeV), P^da/dPT (pb-GeV), for 
I cos^l < 0.77 for the remnant-jet tag without the VDM subtraction. The hatched areas 
are the predictions by the VDM Monte-Carlo simulation. The histograms are predictions 
by the WHIT 1-6 parametrizations. The solid one is WHITl, the dashed one WHIT2, 
the dot-dashed one WHITS, the dotted one WHIT4, the scarce-dotted one WHITS, and 
the short-dashed one WHIT6. Two values of P™"'s were used, i.e., (a) 2.0 and (b) 2.5 
GeV. 
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Figure [1], R. Enomoto et al., Physics Letters B. 
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Figure |2], R. Enomoto et al., Physics Letters B. 
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Figure |3|, R. Enomoto et al., Physics Letters B. 




Figure R. Enomoto et al., Physics Letters B. 
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